Mon. Not. R. Astron. Soc. 000, 000-000 (0000) 



Printed 28 March 2013 



(MN LMftX style file v2.2) 



The spectra of jet bases in FRI radio galaxies: implications 
for particle acceleration 



CO R.A. Laing *\ A.H. Bridle 2 

T ^ i 1 European Southern Observatory, Karl-Schwarzs child- Strafie 2, D-85748 G arching -bei-Milnchen, Germany 
1 2 National Radio Astronomy Observatory, 520 Edgemont Road, Charlottesville, VA 22903-2475, U.S.A. 

<N ■ 



Received 



ABSTRACT 

We present accurate, spatially resolved imaging of radio spectra at the bases of jets 
in eleven low- luminosity (Fanar off- Riley I) radio galaxies, derived from Very Large 
Array (VLA) observations. We pay careful attention to calibration and to the effects 
of random and systematic errors, and we base the flux-density scale on re cent mea- 
surements of VLA primary amplitude calibrators by iPerlev fc Butlerl (|2Q13| ) . We show 
images and profiles of spectral index over the frequency range 1.4 - 8.5 GHz, together 
with values integrated over fiducial regions defined by our relativistic models of the 
jets. We find that the spectral index a, defined in the sense I(y) oc v~ a , decreases with 
distance from the nucleus in all of the jets. The mean spectral indices are 0.66 ± 0.01 
where the jets first brighten abruptly and 0.59 ±0.01 after they recollimate. The mean 
change in spectral index between these locations, which is independent of calibration 
and flux-density scale errors and is therefore more accurately and robustly measured, 
is — 0.067 ± 0.006. Our jet models associate the decrease in spectral index with a bulk 
deceleration of the flow from « 0.8c to <0.5c. We suggest that the decrease is the 
result of a change in the characteristics of ongoing particle acceleration. One possible 
acceleration mechanism is the first-order Fermi process in mildly relativistic shocks: 
in the Bohm limit, the index of the electron energy spectrum, p = 2a + 1, is slightly 
larger than 2 and decreases with velocity upstream of the shock. This possibility is 
consistent with our measurements, but requires shocks throughout the jet volume 
rather than at a few discrete locations. A second possibility is that two acceleration 
mechanisms operate in these jets: one (with p = 2.32) dominant close to the galactic 
nucleus and associated with high flow speeds, another (with p = 2.18) taking over at 
larger distances and slower flow speeds. 

Key words: galaxies: jets - radio continuum:galaxies - acceleration of particles 



1 INTRODUCTION 

Jets from active galactic nuclei (AGN), microquasars and 
gamma-ray bursts can accelerate electrons over enormous 
ranges of energy, generating electr omagnetic radiatio n over 
the entire observable spectrum fe.g. lAbdo et al.ll201lh . They 
may even produce the protons and h eavier nuclei whic h 
make up the highest-energy cosmic rays (|Abreu et al.ll20l"ol ) . 
The mechanisms by which jets accelerate particles to rela- 
tivistic energies are poorly understood, however. The form 
of the electron energy spectrum inferred from observations 
of their synchrotron emission carries information about ac- 
celeration and loss processes. The optically-thin synchrotron 
spectrum observed at radio frequencies typically has an ap- 
proximately power- law form I(y) oc with spectral in- 
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dex a w 0.5 — 0.75. This corresponds to a particle en- 
ergy spectrum n(E)dE = noE~ p dE with p = 2a + 1 « 
2 - 2.5. At low energies where loss processes are unim- 
portant, the slope of this spectrum should be character- 
istic of the acceleration mechanism. For example, first- 
order Ferm i acceleration in strong, non -relativistic shocks 
([Bell Il978l : iBlandford fc Ostrikerl 119781 ) has p = 2 (a = 
0.5), whereas t he asymptotic value of p for ultrarelativisti c 
shocks is 2.23 (|Kirk et alJl200d : iLemoine fc Pelletieill2003h . 
Mildly rela tivistic shocks can produc e a variety of spec- 
tral slopes ([Summerlin fc Baring! [2012h . Other mechanisms 
which have been sugg ested for electron ac c eleration in jets 
inclu d e she ar ([Stawarz fc Ostrowskil 120021 : iRieger fc Duffvl 
12004 l2006|y second-orde r Fermi acceleration by turbu- 
lence (e.g. [Petrosianl l2012h and magnetic reconnection (e.g. 
iBirk fc Le sch 2000). The dominant acceleration mechanism 
may change both locally and globally with position in the 
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jets, for example if the flow decelerates from supersonic 
to subsonic, entrains clumps of denser material, develops 
significant shear or becomes turbulent. Accurate, spatially- 
resolved measurements of the synchrotron spectra from 
AGN jets can, in principle, differentiate between accelera- 
tion mechanisms and show us which processes dominate in 
different regions of the jets. 

The kiloparsec-scale base s of jets in Fanar off- Riley 
Class I radio galaxies (FRI; iFanaroff Rilevl 1 19741 ) are 
known to emit at all frequencies from the radio to the soft 
X-ray regimes. The integrated spectra of the jet bases can 
usually be fit by broken power laws with a ~ 0.5 — 0.7 at low 
frequencies, steepening to a 0.7 — 1.4 in the X-ray band 
(lEvans et alj|2005l: lHardcastle, Birkinsha w Sz Worralllf 
Hardcastle et alj|2005l: lLanz et al.ll201ll; IPerlman ^Wilson 
20051: IPerlman et al l l20ld : IWorrall et alT 



Table 1. Source names, redshifts and linear scales. (1) name as 
used in this paper; (2) alternative names; (3) redshift; (4) linear 
scale in kpcarcsec -1 ; (5) reference for redshift. 



Harwood fc H ardcastle l2012h . 
ization (e.g. I Permian et al.l Il999l 
the emission mechanism at optical 



Ima ging of linear polar- 
20!Ch confirms that 
frequencies is indeed 
synchrotron radiation, as is well established for the radio 
band. The X-ray spectrum is also in accord with this 
interpretation. The corresponding energy distribution 
(also a broken power law) is often referred to as a "single 
electron population" in the sense that it is continuous 
between the presumed, but usually unobservable, lower and 
upper limits. Obse rvations of the c l osest r adio galaxies, 
CenA and M87 (lHardcastle etHI 120071 : IWorrall et all 
120081 : IPerlman fc Wilson! 1 20051 ). show that there are both 
small-scale inhomogeneities and large-scale gradients in the 
broad-band spectra. 

An important implication of the detection of extended 
X-ray emission from AGN jets is that ongoing, distributed 
particle acceleration is necessary to counterbalance syn- 
chrotron cooling. The inferred synchrotron lifetimes of the 
radiating electrons are typically hundreds of years, corre- 
sponding to light-travel distances similar to, or smaller than, 
the jet radii. The acceleration mechanism(s) must therefore 
act over much of the le ngth of a jet, although not necessarily 
over its entire volume (jPerlman &; Wilson! 1 20051 ). Accelera- 
tion cannot be restricted to one or two discrete locations 
such as jet-crossing shocks. 

Precise synchrotron spectral indices at radio frequen- 
cies were de t ermin ed for the jets in two FRI sources by 
lYoung et al.l (|2005h . Combining these results with others 
from the literature, they found (a) = 0.55, apparently in- 
consistent with the value of a = 0.5 (p = 2.0) predicted 
for first-order Fermi acceleration at strong, non-relativistic 
shocks. 

We are studying the kinematics and dynamics of jets in 
FRI radio galaxies by fitting deep, high-resolution VLA ob- 
servations with model brightness dis tributions for symmet- 
rical, deceleratin g relat i vistic flows (lLaing &; Brid le 2002a; 



Canvin fe Laingj 12004); ICanvin et al.1 120051 : lLaing et all 



2006a: lLaingfc Bridle 2012h. In the course of this work, we 



have acquired well calibrated and uniformly imaged data at 
multiple frequencies for the jets in ten FRI radio galaxies. 
These observations, supplemented with additional data from 
the VLA Archive, are ideal for the study of jet base spec- 
tra over the frequency range 1.4 - 8.5 GHz at high spatial 
resolution and we have publish ed spectral-index i m ages for 
many of the sources elsewhere (|Laing et al.ll2006al lbl J2008U 
2011). Data of similar quality, reduced in a similar way, are 



Source 


Alternate 


z 


scale 


Ref 


NGC 193 


UGC408, PKS 0036+03 


0.0147 


0.300 


4 


NGC315 




0.0165 


0.335 


6 


3C31 


NGC 383 


0.0170 


0.346 


5 


0206+35 


UGC1651, 4C 35.03 


0.0377 


0.748 


2 


0326+39 




0.0243 


0.490 


3 


0755+37 


NGC 2484 


0.0428 


0.845 


4 


3C 270 


NGC 4261 


0.0075 


0.154 


6 


M84 


3C 272.1, NGC 4374 


0.0035 


0.073 


6 


3C 296 


NGC 5532 


0.0247 


0.498 


3 


1553+24 




0.0426 


0.841 


1 


3C449 




0.0171 


0.347 


1 



R eferences: 1 Ide Vaucouleurs et"aH (Il99lh: 2lFalco et al.l (119991): 
3 iMiller et al.l (l2002h: 4 | Qgando et al.l (|2QQ8h : 5 ISmith et al.l 
( 2000); 6 Trager et all (l200oh . 



also available for the FRI radio galaxy 3C 449 whose jets are 
too close to the plane of the sky to model using our methods, 
but w hich was studied for different reasons bv lGuidetti et al.l 
(|2010h . 

Our initial work on the jet spectra in 3C3 1, NGC 315 
and 3C 296 ( summ arized in lLaing et al]l2008ah confirmed 
I Young et all (|2005l Vs result that the jet base spectral indices 
can be significantly higher than 0.5 and consistently found 
a value of a « 0.62 near the flaring points where the jets 
brighten abruptly. In a detailed study of the spe ctrum of the 
bright er jet in the giant radio galaxy NGC 315 (|Laing et al.l 
l2006bl ). we found that the jet spectral index decreases with 
distance from the AGN, contrary to the expectations of sim- 
ple synchrotron- loss models, and that there is significant 
transverse spectral structure, with a higher spectral index 
(steeper spectrum) on the jet axis than at its edges. 

In the present paper, our aims are: 

(i) to present a coherent analysis of the jet base spec- 
tra for the whole sample, with carefu l attention to the 
flux-density scale (jPerlev fe Butlerll2013l ). calibration errors, 
noise measurement and correction for contaminating emis- 
sion from surrounding structures; 

(ii) to establish whether the flattening of the spectrum 
with increasing distance from the AGN that we found in 
earlier work is typical of FRI jets; 

(iii) to measure the spectral indices at physically mean- 
ingful locations defined by the geometry, velocity and proper 
emissivity of the jets, to determine the dispersion between 
sources and to search for relations between spectrum and 
flow parameters; 

(iv) to compare these measured values of spectral index 
with the predictions of particle-acceleration models. 

Section [2] presents the sample, gives references to observa- 
tions and basic data reduction (including a discussion of the 
flux-density scale) and describes the additional image pro- 
cessing we have performed to derive the final spectral-index 
images. Section [3] outlines our jet-modelling method and de- 
fines the standard regions for which we tabulate integrated 
spectral indices. In Section 2] we present the total-intensity 
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and spectral-index images, show profiles along the jet axes 
and tabulate integrated values at fiducial locations. We dis- 
cuss the results in Section [5] and summarize our conclusions 
in Section [6] Appendix [A] outlines the data reduction for 
archival observations of two sources which have not been 
presented elsewhere and Appendix iBl gives notes on individ- 
ual sources, primarily to document instrumental issues and 
to reference earlier work on spectral imaging. 



2 OBSERVATIONS AND DATA REDUCTION 

2.1 The sources 

Our modelling technique requires us to be able to image the 
jets on both sides of the nucleus with good transverse resolu- 
tion at high signal-to-noise ratio, and for their intensity ratio 
close to the nucleus to be >5:1 (Section[3]). This led to the 
selection of ten sources. We have added 3C 449, which sat- 
isfies the first criterion, but not the second. The sources are 
listed in Table [T] together with their redshifts and the cor- 
responding linear scales for a concordance cosmology with 
Q M 0.3, Q A 0.7 and H 70 krns" 1 Mpc -1 . 

2.2 VLA images 

2.2.1 Observations and references 

The observations were made using various combinations of 
frequencies and array configurations in the 1.3 - 1.7 GHz, 
4.5 - 5.0 GHz and 8.3 - 8.5 GHz bands of the VLA. We 
have published comprehensive descriptions of the observa- 
tions and data reduction for almost all of the images used in 
this paper. The only exceptions are the lower- frequency im- 
ages for two sources, 0326+39 and 1553+24. These have not 
been published elsewhere and the reduction of VLA archive 
observations is described in Appendix [A] All of the images 
used in the present paper are listed in Table [2] with refer- 
ences to descriptions of the observations and reductions. We 
summarize some general points below. 

2.2.2 Calibration and imaging 

The data were calibrated, imaged and self- calibrated using 
standard techniques in the MPS package. 3C 286 was used 
as the primary amplitude calibrator if it was observable 
at sufficiently high elevation during an observation, other- 
wise 3C48 was used. We discuss the flux-density scale in 
Section 12.2.31 below. Our observations were designed more 
for accurate imaging and polarimetry than for flux-density 
measurement, and we did not determine gain- elevation or 
opacity corrections independently. The transfer of the flux- 
density scale was always done with secondary and primary 
calibrators observed at similar elevations, however, so we do 
not expect large errors from these effects over our frequency 
range. 

All images were made using data from multiple con- 
figurations of the VLA. We did not attempt to match the 
spatial-frequency coverage of the datasets in detail, but 
rather to ensure that the relevant range of spatial scales 
was covered fully. All of the visibility datasets for the indi- 
vidual frequencies were imaged with scaled maximum base- 
lines, Gaussian tapers and data weighting to give the same 



resolution and were deconvolved and restored with the same 
Gaussian beam. We used three different deconvolution algo- 
rithms (see the references in Table [2}, depending on the 
complexity of the source structure: conven tional (single- 
resolution) clean; multi-res olution CLEAN (iGreisen et al.l 
2009) and maxim um entropy (ICornwell &; Evansll 19851 , used 
as described by iLeahv &; Perlevl Il99lh . Data were usually 
taken in two frequency channels simultaneously. At fre- 
quencies above 4.5 GHz, the two channels were imaged to- 
gether if they were contiguous. The channels in the 1.3 - 
1.7 GHz band (typically more widely spaced in frequency) 
were imaged separately. In a few of these cases where the 
centre frequencies of two channels were less than 100 MHz 
apart, the images were averaged after deconvolution. 

The registration of the images was assured during the 
imaging process by shifting the visibility data so that the 
unresolved central component (or "core" ) was always aligned 
on the centre of a given pixel (this was in any case required 
for accurate deconvolution). 



2.2.3 Flux-density scale 

The accuracy and repeatability of the relative flux-density 
scale across our frequency range are clearly important for 
this study. We initially used the values of flux density and 
visibility models for 3C 28 6 and 3C48 provided with the 
AlPS software distribution. IPerlev fe Butler] (|2Q13h recently 
carried out a comprehensive investigation of the VLA flux- 
density scale. They showed that the flux density of 3C 286 
was stable over the period of our observations, but that of 
3C 48 varied significantly. They derived improved flux den- 
sities for bot h calibrators on th e absolute scale established 
by WMAP (jWeiland et alJfeOllh . extrapolated to lower fre- 
quencies using a thermophysical model of Mar s. They also 
i ncorp orated flux-density measurements from iBaars et al.l 
(|l977l ) at frequencies <3 GHz. We rescaled all of the images 
used in the present paper, multiplying by th e ratio of the flux 
density of the primary calibrator given bv lPerlevfc Butlerl 
to that assumed in the initial data reduction. Specif- 
ically^ ^we_as^d__tl^ expressions in Tables 10 and 
11 of Per lev fe Butlerl {2013) for our observing frequencies, 
interpolating linearly to the relevant epoch for 3C 48. Data 
taken over periods of a few months to several years were 
combined to derive the final images, so the variability of 
3C 48 between individual observations was a potential con- 
cern. For the relevant combinations, its maximum variability 
amplitude was <1%, however, so we ignored this complica- 
tion. 

The average effect of the recalibration determined by 
IPerlev fe Butlerl {2013) is to increase the measured flux den- 
sities at 1.3 - 1.7 GHz by ^2% and to decrease those at 4.5 - 
5.0 GHz and 8.3 - 8.5 GHz by similar amounts. The result- 
ing change in spectral index is between +0.03 and +0.04. 
Sources are differently affected by the variability of 3C48 
and (very slightly) by inconsistencies in the flux densities 
assumed for initial calibration. The individual correction fac- 
tors are listed in Table [2] 



1 The exceptions are early observations or those taken during 
phased-array VLBI runs. 
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Table 2. Parameters of the VLA images used to derive spectral indices. (1) source name; (2) resolution (FWHM), in arcsec; (3) rotation 
of image with respect to the sky, in deg; (4) centre frequency, in MHz; (5) VLA configurations used to make the im age (H denotes a 
hybrid configuration); (6) correction factor applied to published observations to bring them onto the flux-density scale of lPerlev fc Butlen 
(|2013l ); (7) off-source noise level, in /xJybeam -1 ; (8) noise level close to the core, in /iJybeam -1 ; (9) reference for image; (10) method 
used to derive spectral-index images and averaged values (2 for two- frequency; PL for a power-law fit). Cols 11-13 give the parameters 
for lobe subtraction, in arcsec. (11) minimum width of reference region; (12) maximum width; (13) smoothing kernel. Cols 11 - 13 are 
blank if no lobe subtraction was done. 



Source FWHM 
arcsec 

NGC193 1.6 



NGC315 5.5 



3C31 1.5 



0206+35 1.2 



0326+39 1.5 



0755+37 1.3 



3C270 1.65 



M84 1.65 



3C296 1.5 



1553+24 1.5 



3C449 1.25 



PA 



175.2 



-41.5 



-70.3 



-41.0 



3.0 



158.5 



3.5 



3.0 



-125.9 



-46.0 



80.4 



Freq 
(MHz) 



Conf 



Corr 



^off Ccore 

MJyb- 1 



Reference 



Fit Min Max Smooth 
arcsec 



4860.1 


ABCD 


0.991 


7.5 


15 


1365.0 


ABC 


1.029 


36 


72 


4985.1 


ABCD 


0.955 


15 


41 


1665.0 


BC 


1.019 


37 


120 


1485.0 


BC 


1.023 


37 


37 


1413.0 


ABC 


1.026 


38 


46 


1365.0 


BC 


1.026 


41 


41 


8440.0 


ABCD 


0.956 


11 


18 


4985.0 


BCD 


0.964 


12 


22 


1636.0 


ABCD 


1.013 


72 


150 


1485.0 


ABCD 


1.018 


73 


150 


1435.0 


ABCD 


1.020 


78 


110 


1365.0 


ABCD 


1.023 


76 


130 


4860.1 


BC 


0.988 


12 


120 


1425.0 


AB 


1.026 


19 


210 


8460.1 


ABCD 


0.967 


6.9 


14 


4885.1 


HCD 


0.980 


22 


51 


1425.0 


ABCD 


1.025 


44 


99 


4860.1 


BCD 


0.980 


7.8 


13 


1425.0 


ABC 


1.022 


20 


36 


4860.1 


ABCD 


0.982 


21 


37 


1388.5 


ABC 


1.022 


21 


50 


4860.1 


ABC 


0.980 


15 


39 


1413.0 


AB 


1.022 


140 


150 


8460.1 


ABCD 


0.971 


9.5 


9.5 


1406.75 


ABCD 


1.021 


27 


48 


8460.1 


ABC 


0.974 


8.5 


13 


4860.1 


BCD 


0.980 


22 


30 


1385.1 


ABC 


1.023 


35 


42 


8385.1 


ABCD 


0.987 


14 


24 


4985.1 


ABCD 


0.982 


18 


28 


4685.1 


ABCD 


0.983 


17 


27 


1485.0 


ABCD 


1.021 


35 


45 


1465.0 


ABCD 


1.021 


48 


58 


1445.0 


ABCD 


1.021 


21 


31 


1365.0 


ABCD 


1.023 


37 


47 



Laing et al. (2011) 2 
Laing et al. (2006b) PL 



12 18 



Laing et al. (2008b) PL 



Laing et al. (2011) 



Canvin & Laing (2004) PL 
Appendix [Al 



Laing et al. (2011) 



10 



30 45 



Laing, Guidetti & 2 14.4 18 

Bridle, in preparation 



Guidetti et al. (2010) PL 
Feretti et al. (1999) 



1.5 



1.0 



3.0 



2.0 



1.5 



Laing et al. (2011) 2 6 9 

Laing et al. (2006a) 2 17.7 19.2 1.5 



Canvin fe Laing (2004) PL 
Appendix [Al 



2.3 Lobe subtraction 

Seven of the sources in our sample have lobes surround- 
ing their jets on scales of interest to this work (NGC 193, 
0206+35, 0326+39, 0755+37, 3C270, M84 and 3C296). In 
six of these, the surface brightness of the lobes is significant 
at arcsecond resolution (the exception is 0326+39). In order 



to isolate and measure the spectrum of the jet emission, we 
need to remove the contribution of the lobe, which typically 
has a steeper spectru m than the jets. To do this, we used the 
method described by lLaing fe Bridle! (|2012l ). assuming that 
the lobe intensity varies smoothly and slowly across the jet. 
We first located the edges of the flatter-spectrum jet emis- 
sion on images of spectral index and/or by using spectral 
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tomography (Kat z-Stone fc RudnickHl997T ). We then defined 
two background regions parallel to the jet axis and just out- 
side the maximum extent of the jet emission, smoothed the 
background brightness distributions parallel to the jet axis 
with a boxcar function and interpolated linearly between 
them under the jet. The locations of the background regions 
and the widths of the smoothing kernels are given in Table [2] 



2.4 Error model 

Errors in measured brightness can be divided into three 
types: random errors which vary across individual images; 
calibration errors which multiply all pixels within a single 
image by the same factor but are uncorrelated between im- 
ages and scale errors which affect all observations at given 
frequency in the same way. 

We took the random component to have additive and 
multiplicative sub-components. For the additive part, we as- 
sumed that errors in total intensity have a Gaussian distri- 
bution with zero mean and rms 07 in the image plane and 
that they are independent on scales larger than the Gaussian 
restoring beam, which has an area of 7r/ 2 (21n2)- 1 , where / 
is the FWHM. When computing spectral-index images and 
determining blanking levels, we took 07 to be the off-source 
noise level cr Q ff averaged over several areas well clear of any 
emission as given in Table [2] All of the sources in this study 
have bright, unresolved cores, so the effect of residual cali- 
bration errors is to increase the noise level at the jet bases. 
When computing spectral indices from flux densities inte- 
grated over different regions of the jets we therefore followed 
the following recipe. 

(i) We used no data from points closer than 2f to the 
core. 

(ii) We only included pixels in the summation if the ran- 
dom error on the spectral index at that point determined 
using the off-source noise levels cr ff was <0.03 (this auto- 
matically ensured that we integrated over the same points 
at all frequencies). 

(hi) For points farther than 4/ from the core, we took 

o~i — CT ff. 

(iv) We adopted a different noise level a core for pixels lo- 
cated between 2f and 4/ from the core. cr CO re, which is also 
given in Table [2] was estimated from rectangular regions 
displaced perpendicular to the jet axis on either side of the 
core. 

We also included a random error contribution which scales 
with surface brightness, as might be expected for some types 
of calibration or deconvolution artefacts. Such effects are 
potentially quite heterogeneous, ranging from offsets in the 
zero level due to missing short-spacing coverage which af- 
fect an entire image, through coherent ripples with various 
wavelengths to spotty structures caused by the failures of 
the clean algorithm. For our images, they become compa- 
rable with the additive component only after averaging over 
many beamwidths. On the basis of a comparison of images at 
the same or neighbouring frequencies, we therefore crudely 
modelled this type of error as a 1% rms multiplicative effect 
applied to integrated flux densities only (i.e. assumed to be 
coherent over the integration region). We also carefully re- 
viewed all of the total intensity and spectral-index images 



and identified only one source, M84, where coherent cali- 
bration/deconvolution errors on scales larger than the inte- 
gration regions must be significant. 0206+35 shows artefacts 
at a lower level, consistent with our error prescription. Both 
cases are discussed in Appendix 151 

We next considered multiplicative calibration offsets. 
Errors in the gain calibration process may be partially cor- 
related between datasets for the same source, given that 
observations at different frequencies were often interleaved. 
Nevertheless, all of the images were made from at least two 
and often many more individual datasets, and we believe 
that errors in the transfer of the flux-density scale from pri- 
mary calibrator to target are essentially independent. Addi- 
tional errors may be introduced during combination of data 
from different times and/or array configurations. Although 
we always forced the mean modulus of the complex gain 
correction to be unity during self-calibration, the simulta- 
neous adjustment of the relative amplitude scale of the two 
datasets and the flux density of the compact core, which ma y 
vary with time (|Laing &; Bridle! l2002al : lLaing et al.l l2006al ) , 
was not always fully constrained at the higher frequencies, 
and additional errors were introduced. Our best estimate, 
derived from internal consistency checks and in particular 
from a comparison of independent measurements with the 
same array configuration (s) at similar or identical frequen- 
cies, is that the appropriate rms calibration error for an 
individual image is 3% at all of our frequencies, or more for- 
mally that the measured brightness is kl, where log/c has 
a Gaussian distribution with rms 0.0128. We discuss the ef- 
fects of inaccuracies in this estimate at appropriate points 
below. We treat this type of error as uncorrelated between 
sources. 

Finally, errors in the flux densities of the primary am- 
plitude calibrators introduce systematic scale errors for all 
observations at a given frequency. The absolute accuracy 
of the new VLA flux-density scale is thought to be better 
than 2.5% r ms at 1.3 - 1.7 GHz a nd 1% rms at our higher 
frequencies (|Perlev fe Butler! I2013T ). We are concerned with 
the relative accuracy over our frequency range, which must 
be somewhat better than this. 

2.5 Spectral-index fitting 

If images were available at only two frequencies v\ and 
1/2, we calculated the spectral indices directly as ot v v \ — 
ln[/(z/i)//(i/2)]/ln(z/2/^i) and determined the random error 
cr r from the rms errors in I (following the prescription in 
Section 12. 4p by standard error propagation. The rms cali- 
bration error cr c then just depends on the ratio since 
we assume the same fractional error at all frequencies: 

_ 0.03 x 2 1/2 
ln(z/ 2 /z/i) 

Deviations from power-law spectra are of physical in- 
terest as indicators of synchrotron losses and potentially as 
diagnostics of the acceleration mechanism. For four sources 
(3C31, 0326+39, 1553+24 and 3C449), we have data in all 
of the 1.3 - 1.7, 4.5 - 5.0 and 8.3 - 8.5 GHz bands and could 
potentially have addressed this issue. For each source, we 
chose one frequency per band (^w, ^mid, ^high) and plotted 
a^f d h against a£™j^ (colour- colour plots; iKatz-Stone et al.l 
1993) for the binned profiles from Section H~2l (see below). 
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For all four sources, we found that the points were tightly 
clustered around single locations in the colour-colour dia- 
gram, consistent with the effect of a calibration error. From 
our error model (Section 12. 4[) , we predict an rms for the 
spectral- index difference ctf^f* — a^jf of 0.10; we found 
an rms of 0.12 and a mean of +0.05 for the four sources, 
with one (3C31) showing spectral flattening with increasing 
frequency. We conclude that any apparent spectral curva- 
ture observed over our relatively narrow frequency range 
is entirel y consistent with t he effects of calibration errors. 
Fig. 13 of lLaing et al.ll2008bl shows integrated spectra which 
reinforce this point. 

Given our lack of evidence for real spectral curvature, 
we chose to fit single power laws to spectra with three or 
more frequencies and determined the spectral indices by 
least-squares fitting, using weights derived from the random 
and systematic errors, summed in quadrature. We then took 
a c to be the error in the limit of zero random noise and 
cr r = (cr| t — <7c) 1//2 , where am is the error estimate from the 
fit. 

For consistency, all of the displayed spectral-index im- 
ages (whether determined from two frequencies or more) are 
blanked wherever a r > 0.03, assuming constant noise levels 
o~i = cr ff across the maps to calculate cr r (i.e. without mul- 
tiplicative random or systematic errors). This criterion is 
also part of that used to select the points included in the 
integrated flux densities (Section 12.41 above). 



Table 3. Fiducial distances derived from jet models, projected 
on the plane of the sky, and details of jet detections at higher fre- 
quencies. (1) name; (2) recollimation distance zo] (3) start of high- 
emissivity region (flaring point), z e \; (4) end of high-emissivity 
region z e o] (5) start of deceleration, z v i; (6) end of deceleration 
Zvo; (7) maximum length of detected X-ray or optical emission 
zox for the main jet; (8) notes and references for high-frequency 
detections. All distances are in arcsec. Only the recollimation dis- 
tance and the position of the flaring point are given for 3C449: 
these are the two parameters that can be determined directly 
from the brightness distribution without a model. 



Name Distances / arcsec Note 







z e i 


ZeO 


Zvl 


Z V 






NGC193 


22.2 


3.3 


7.0 


0.0 


13.9 


5 


1 


NGC315 


79.5 


5.8 


17.7 


15.1 


36.8 


30 


2 


3C31 


8.5 


2.8 


8.5 


5.4 


9.5 


8 


3 


0206+35 


4.5 


0.7 


1.8 


1.5 


3.5 


2 


4 


0326+39 


18.2 


1.8 


5.1 


3.0 


12.4 






0755+37 


9.4 


1.0 


9.0 


2.5 


18.7 


4 


5 


3C 270 


34.0 


9.0 


25.3 


13.0 


28.7 


32 


6 


M84 


17.0 


2.8 


6.8 


2.7 


13.2 


4 


7 


3C 296 


27.5 


3.0 


13.0 


8.9 


9.8 


10 


8 


1553+24 


6.5 


0.5 


1.7 


3.7 


4.2 


1 


9 


3C 449 


11.0 


4.5 













3 JET MODELS 

Our models assume that the j ets are symmetrical , 
axisymmetric and re l ativistic (ILaing &; Bridle! 2002a ; 



Canvin Laingl [20041 : 
2006al: ILaing fc Bridl. 



, Canvin et 
2012). 



Lams & Bi 
alj [2005; 



Laing et all 



Relativistic aberration not 
only causes the approaching (main) jet to appear brighter 
than the receding counter-jet, but also leads to differences 
in observed linear polarization, as the angles to the line of 
sight in the rest frame of the emission are different on the 
two sides. Simultaneous fitting of Stokes /, Q and U then 
allows us to break the degeneracy between velocity and in- 
clination. The basic steps in our method are: 

(i) construct a parameterized model of geometry, velocity, 
emissivity and field-ordering in the rest frame; 

(ii) calculate the observed-frame emission in /, Q and U ; 

(iii) integrate along the line of sight and convolve with 
the observing beam; 

(iv) evaluate % 2 between model and observed images and 

(v) optimize the model parameters. 

In order to make physically meaningful comparisons be- 
tween spectral indices in different jets, we defined integration 
regions based on three aspects of our model fits: geometry, 
proper emissivity and velocity. In Fig. [1] we illustrate these 
regions using the observed main-jet brightness distribution 
and model fits for NGC 315 (Laing & Bridle, in preparation). 
The relevant aspects of our model fits are as follows. 

Geometry We divide a jet into two regions based on the 
geometry of its outer isophotes: a flaring region, where the 
flow first expands with increasing opening angle and then 
recollimates and a conical outer region (Figs [TJi and b). 
We define zq to be the distance between the AGN and the 



References. Unless noted explicitly, these describe X-ray observa- 
tion^ 

ll Kharb et~aD (120121) 

2 Worrall et al l (l2Q07h 

3 Hardcastle et al .1 (|2002h : optical and mid-IR emission ha ve also 
been detected from the main jet (jLanz et al. 2011; Cr oston et al.l 
l20Q3h 

4 lWorrall. Birkinshaw Hardcastle] (|200lh 

5 IWorrall et al.l (|200ll ~ the inner 2.6 ar csec of th e main jet has 
b een detected a t optic al wavelengths (|Parma et al.ll2003h 
6 lWorrall et al.l (|2010T ): the inner 20 arcsec of the counter-jet was 
d etected in the same observation. 

7 [Harris et al.l (|2002|) 

8 lHardcastle et all §2005); optical emission from the main jet is 
a lso reported in this reference . 

9 lParma etHI (2003); this is a detection of optical emission. 



boundary between the regions, measured along the axis and 
projected on the plane of the sky. 

Proper emissivity We model the quantity noB 1+a , 
where no is the normalizing constant in the electron energy 
distribution and B is the rms total magnetic field, both mea- 
sured in the rest frame, no and B cannot be separated using 
observations of synchrotron emission without additional as - 
sumptions, for example flux-freezing (jLaing fc Bridldl2004h 
or equipartition between particle and field energy. We refer 
to noB 1+a , loosely, as "the emissivity"; in fact it is multi- 
plied by a constant and by factors dependent on field geom- 
etry (which we ignore in the pr esent discu ssion) to give the 
true emissivities in /, Q and U (|Laing|l2002l ). The emissivity 
profile along the jet axis is divided into 3 or 4 sections, in 
each of which the emissivity has a power-law dependence on 
distance from the nucleus. All of the jets have a well-defined 
flaring point marked by an abrupt increase in brightness. We 
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model this as a step in emissivity and a change in power- 
law slope. The flaring point marks the start of the high- 
emissivity region, which is characterized by bright, knotty, 
non-axisymmetric substructure and which is usually obvious 
by eye in the brightness distribution of the main jet. The end 
of this region is again marked by a step in emissivity and/or 
a change (usually flattening) of slope. These features are il- 
lustrated for NGC315 in Figs [TJc) and (d). The positions 
of the start and end of the high-emissivity region from the 
AGN, measured from the nucleus along the projection of the 
jet axis on the sky, are z e i and z e o, respectively. 

Velocity The on-axis velocity is modelled as initially con- 
stant, decreasing linearly with distance in the deceleration 
region and then varying, again linearly but much more 
slowly, at large distances (we allow either acceleration or de- 
celeration in the fits) . Deceleration usually begins within the 
high-emissivity region and ends before recollimation. The 
start and end of the deceleration region, again measured 
from the AGN along the projected jet axis, are z w \ and z v o, 
respectively. 

A full tabulation of fitted model parameters will be given 
by Laing & Bridle (in preparation). In the present paper, 
we use only the fiducial distances and (in Section f 5 . 2 1) , the 
velocity profiles. We list the angular distances zq, z e i, z e o, 
z v ± and z v o in Table [3] These quantities define our integra- 
tion regions. Other aspects of our models, such as the shapes 
of the streamlines and boundary surfaces between regions, 
the transverse variations of velocity and emissivity and the 
distribution of magnetic-field component ratios are less rel- 
evant to the present discussion. 

We defined four integrated spectral indices, as follows. 

High-emissivity, ahe We integrated over the part of the 
high-emissivity region that is farther than 2f from the core 
(Figs [It and d), stopping the integration at the end of the 
deceleration region if the high-emissivity region extends be- 
yond it. 

Deceleration, c^decei We summed over the deceleration 
region, stopping at recollimation if this occurs before the 
end of deceleration. 

Coasting, a CO ast We integrated from the end of the de- 
celeration region for 4/ or until recollimation, whichever is 
closer. 

Recollimation, a r ecoii The integration started at recolli- 
mation and extended for 4f into the outer region. 

The high-emissivity and deceleration regions usually over- 
lap, but the definitions were designed to ensure that the 
coasting and recollimation regions are disjoint with each 
other and with the remaining two regions. The resolution 
of our observations is too low to measure reliably the spec- 
tral index of emission any closer to the AGN except in the 
main jet of 3C 270 (Appendix [5]). 

We integrated the flux densities over all pixels satisfying 
the criteria of Section 12.41 requiring in addition that the to- 
tal number of unblanked pixels in the region, N > 2A^beam, 
where A^beam = 7r(//Ax) 2 (2 In 2) _1 is the beam area in pix- 
els of size Ax. We then calculated spectral indices and errors 
as in Section [231 




Distance from AGN / arcsec 



Figure 1. Jet models and fiducial regions for spectral indices, 
(a) The main jet of NGC315 at 2.35-arcsec FWHM resolution, 
indicating the area over which a reC oll was evaluated, (b) Sketch 
showing the flaring and outer regions of the jets in NGC 315 to 
the same scale as panel (a), but with no bend in the jet. (c) Lon- 
gitudinal profiles of rest-frame emissivity derived from our model 
of NGC 315 and plotted against distance from the AGN projected 
on the plane of the sky. The high-emissivity region is marked, (d) 
Grey-scale of the main jet at a resolution of 0.4 arcsec FWHM, 
showing the high-emissivity, deceleration and coasting regions, 
(e) Longitudinal profiles of velocity derived from our model of 
NGC 315, with the deceleration and coasting regions marked. In 
panels (c) and (e), the full and dashed lines apply to the axis and 
edge of the jet, respectively. 



4 SPECTRAL-INDEX IMAGES AND 
PROFILES 

4.1 Images 

In Fig. [2] we show false-colour images of spectral index over 
the range 0.4 < a < 0.8 for all of the sources, together with 
grey-scales of total intensity at the same resolution. The 
images have been rotated by the angles given in Table [2] so 
as to have the main jets on the right. 

Except for the unresolved cores, which are partially 
optically thick with a<0.3, the emission typically has 
0.7>a>0.5. Significantly higher or lower spectral indices 
are only found close to the edges of the jets (where small er- 
rors in deconvolution, lobe subtraction or image registration 
can cause large changes in a) or in regions of low signal-to- 
noise ratio. The mean spectral indices of the jets vary over a 
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(a) NGC 193 



(b) NGC 315 



(c) 3C 31 




(j) 1553+24 


(k) 3C 449 







0.4 0.5 0.6 0.7 0.8 



Figure 2. False-colour images of spectral index, a (upper panels) and grey-scales of total intensity (lower panels) at the same resolution. 
a is in the range 0.4 - 0.8 (as indicated by the labelled wedge) and is plotted wherever cr r < 0.03. The images have been rotated by the 
angles given in Table [2] so that the approaching jets are always on the right. On the spectral-index plots, the horizontal lines represent 
distances of 10 arcsec and the restoring beams are shown in the bottom-left-hand corners. Lobe emission has not been subtracted from 
the / images. 
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Figure 3. Profiles of spectral index, a, along the jets. The brighter (approaching) jet is always on the right. The profiles are derived 
from the images shown in Fig. [2] along horizontal lines passing through the nuclei (a different position angle was used for the counter-jet 
of M84, which is misaligned with the main jet). Profiles are only plotted where cr r < 0.03 and the distance from the core exceeds 2 x 
FWHM. The beamwidths (FWHM) are indicated by the horizontal lines centred on distance and a = 0.45. The vertical dashed black 
lines represent the extent of the high-emissivity region and the dotted vertical green lines the start and end of rapid deceleration. The 
vertical full cyan lines show where the jets recollimate. The horizontal red lines plotted at a = 0.5 or a = 0.7 show the maximum extent 
of any detected soft X-ray or optical emission from the jets (Table EJ). 
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Figure 4. Profiles of integrated spectral index, a, along the jets. The spectral indices were calculated from total intensities integrated 
over points with cr r < 0.03 in boxes of length 2FWHM along the jet axis (horizontal in Fig. [2] starting at 2FWHM from the core. Data 
are only plotted if there are more than two independent beam areas above the blanking level in the box. Blue full and red dashed symbols 
represent main and counter-jets, respectively. The fiducial distances are indicated as in Fig. [3] 
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range of ^0.1, with NGC315 showing the flattest spectrum 
and 3C270 the steepest. The spectral indices of main and 
counter-jets are usually very similar, with M 84 and 1553+24 
the only conspicuous exceptions. 

A slight, but significant tendency for the spectral index 
to decrease away from the AGN is apparent from the images 
in Fig. [2] The effect is subtle (Aa<0.1) but consistent. 

4.2 Longitudinal profiles 

In order to investigate the spectral gradients in more de- 
tail, we have plotted profiles of spectral index along the jets 
in two representations: as unaveraged slices along the jet 
ridge-lines to show the level of variation on small as well as 
large scales (Fig. [3]) and as binned averages to reduce ran- 
dom noise and to compare values in the main and counter- 
jets (Fig. [J). The unbinned profiles are plotted for distances 
>2f from the AGN wherever the rms random error on the 
spectral index determined from the off-source noise cr ff (Ta- 
ble [2]) is cr r < 0.03. The binned profiles follow the prescrip- 
tion given in Section [2. 41 

The clear tendency for the spectra to flatten with dis- 
tance from the AGN is confirmed, and the binned profiles 
(Fig. 3|) show in addition that the jet and counter-jet spec- 
tral indices at the same distance from the AGN are usu- 
ally very similar. Of the two exceptions noted earlier, M 84 
(Fig. shows an offset between the spectral indices of its 
jets which results from an instrumental effect (Appendix [Bj 
and 1553+24 (Fig. [4]) has a weak counter-jet with measur- 
able spectral index only over a limited range of distance. 

The fiducial distances from our jet models (projected 
on the sky) are plotted in Figs [3] and [H Most of the spectral 
flattening with distance from the AGN occurs before the jets 
recollimate, and appears to be associated with rapid decelera- 
tion. The clearest examples are the well-resolved and bright 
sources NGC315, 3C31, 3C 270 and 3C296 (Figs Bt>, c, g 
and i). 

High-frequency (>10 14 Hz) emission is seen in many of 
the main jets and in the counter-jet of 3C 270. Although the 
data are heterogeneous, with a mixture of optical and X-ray 
detections and various exposure times, the comparison in 
Fig. [3] shows that the detected high-frequency emission typ- 
ically extends roughly as far as the end of the high-emissivity 
region. 

4.3 Spectral indices at fiducial locations 

Individual integrated spectral indices over the four fidu- 
cial areas are given in Table [4j along with rms random 
errors a r estimated using the prescription in Section 12.41 
The rms calibration errors cr c , which are common to all of 
the measurements for an individual jet, are also tabulated. 
cr r w 0.01 — 0.02 in regions of high signal-to-noise, with sig- 
nificantly larger values primarily in those counter-jets whose 
high-emissivity regions are both faint and close to a bright 
core. 

Histograms of the distributions are plotted in Fig. [5] 
These display the progressive decrease of spectral index 
from high-emissivity through deceleration to coasting re- 
gions. There is no evidence of further systematic spectral 
flattening after recollimation: the distributions for the coast- 
ing and recollimation regions are statistically the same (the 
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Figure 5. Histograms of spectral indices at different locations, as 
given in Table |4] (a) high-emissivity; (b) deceleration; (c) coast- 
ing; (d) recollimation. The shaded areas represent the main jets. 
The vertical dotted lines indicate the unweighted means of the 
distributions, from Table [5] 

significance level of the difference is only 28% according to 
the Kolmogorov-Smirnov test), consistent with very similar 
spectral-index and error distributions. The histogram for the 
high-emissivity region (Fig. [5Jl) is slightly broader than the 
other three; this is expected from the larger random errors 
in the counter-jet spectral indices, and the distribution for 
the main jets alone is narrower. 

Statistics of the distributions for the main and counter- 
jets separately and for the two combined are listed in Ta- 
ble [5] Here, and wherever we quote the error on a mean 
spectral index, we give the standard error derived from 
the dispersion in the measurements. The unweighted means 
are very close to those weighted by the combined (random 
and calibration) errors. The unweighted mean spectral in- 
dex for the high-emissivity region, («he) = 0.662 ± 0.011, 
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Table 4. Integrated spectral indices for high-emissivity, deceleration, coasting and recollimation regions. (1) source name; (2) rms 
calibration error for spectral indices cr c . The remaining columns give the mean spectral indices, (a) and their estimated random errors 
for the four fiducial regions in each jet. Values are not given if one or more of the regions: has too little data of adequate signal-to-noise 
(NGC315; 3C31, 0755+37), is too close to the core (0206+35 and 1553+24) or is undefined (3C449). 



Source cr c High-emissivity Deceleration Coasting Recollimation 

jet cj jet cj jet cj jet cj 

(a) ± (a) ± (a) ± (a) ± (a) ± (a) ± (a) ± (a) ± 



NGC 193 


0.036 


0.64 


0.01 


0.75 


NGC 315 


0.027 


0.65 


0.01 


0.76 


3C 31 


0.017 


0.64 


0.01 


0.64 


0206+35 


0.036 








0326+39 


0.023 


0.67 


0.01 


0.60 


0755+37 


0.036 


0.58 


0.01 




3C 270 


0.036 


0.67 


0.02 


0.71 


M 84 


0.036 


0.64 


0.01 


0.67 


3C 296 


0.034 


0.64 


0.01 


0.62 


1553+24 


0.023 








3C 449 


0.016 


0.68 


0.01 


0.68 



0.05 


0.62 


0.01 


0.63 


0.02 


0.59 


0.01 


0.04 


0.60 


0.01 


0.57 


0.01 


0.55 


0.01 


0.01 


0.63 


0.01 


0.63 


0.01 








0.63 


0.02 


0.56 


0.05 


0.59 


0.03 


0.09 


0.65 


0.01 


0.62 


0.01 


0.60 


0.01 




0.57 


0.01 










0.03 


0.68 


0.02 


0.67 


0.02 


0.64 


0.02 


0.02 


0.61 


0.01 


0.63 


0.01 


0.55 


0.01 


0.01 


0.63 


0.01 


0.62 


0.01 


0.60 


0.01 




0.57 


0.02 






0.60 


0.01 



0.03 



0.59 


0.02 


0.58 


0.02 


0.65 


0.03 


0.55 


0.01 


0.53 


0.01 










0.60 


0.01 


0.60 


0.01 


0.57 


0.04 


0.56 


0.01 


0.56 


0.01 


0.60 


0.01 


0.59 


0.01 


0.57 


0.01 


0.63 


0.02 


0.63 


0.02 


0.64 


0.02 


0.61 


0.01 


0.55 


0.01 


0.59 


0.01 


0.59 


0.01 


0.58 


0.01 


0.57 


0.01 


0.61 


0.04 


0.58 


0.01 


0.67 


0.02 






0.61 


0.01 


0.61 


0.01 



Table 5. Mean and rms spectral indices for the four fiducial re- 
gions, for main and counter-jets separately and for the combina- 
tion. (1) jet class; (2) location; (3) number of jets; (4) unweighted 
mean spectral index; (5) measured spectral index rms; (6) rms 
predicted from the combination of random and calibration errors; 
(7) standard error on the unweighted mean; (8) mean weighted 
by random and calibration errors. 



Jet Region TV Unweighted Wtd 

(«> ^obs ^pred ± (<*) 



Table 6. Differences between spectral indices integrated over the 
standard regions. (1) spectral-index combination; (2) number of 
measurements (main and counter-jets); (3) mean spectral-index 
difference; (4) rms; (5) standard error on the mean; (6) signif- 
icance level of correlation between spectral indices at the two 
locations, from the Spearman rank test. 



Location 


N 


(Aa) 


rms 


± 


sig 


«he - «decel 


15 


0.035 


0.053 


0.014 


77.0 


a he — «coast 


12 


0.077 


0.058 


0.017 


18.0 


^he ^recoil 


15 


0.067 


0.024 


0.006 


99.9 


C^decel C^coast 


15 


0.029 


0.022 


0.006 


>99.9 


C^decel ^recoil 


16 


0.037 


0.027 


0.007 


99.7 


CKcoast — CWoll 


15 


0.005 


0.028 


0.007 


99.9 



Table 7. Mean differences between the spectral indices of the 
main and counter-jets at the four fiducial locations. (1) location; 
(2) mean (a c j — ctj); (3) standard error on the mean; (4) number 
of sources; (5) significance level of correlation between cn c j and 
ckj, according to the Spearman rank test. 



Location (a c j — aj) ± N sig 



Main 


high- 


9 


0.645 


0.028 


0.033 


0.009 


0.650 


CJ 


emiss- 


8 


0.680 


0.054 


0.050 


0.019 


0.668 


Both 


ivity 


17 


0.662 


0.046 


0.042 


0.011 


0.656 


Main 


decel- 


10 


0.619 


0.031 


0.034 


0.010 


0.621 


CJ 


eration 


8 


0.617 


0.035 


0.038 


0.012 


0.618 


Both 




18 


0.618 


0.033 


0.036 


0.008 


0.620 


Main 


coast- 


8 


0.590 


0.026 


0.036 


0.009 


0.590 


CJ 


ing 


8 


0.593 


0.026 


0.040 


0.009 


0.590 


Both 




16 


0.592 


0.026 


0.038 


0.007 


0.590 


Main 


recoll- 


10 


0.582 


0.027 


0.032 


0.009 


0.589 


CJ 


imation 


9 


0.607 


0.036 


0.034 


0.012 


0.605 


Both 




19 


0.594 


0.034 


0.033 


0.008 


0.596 



is significantly steeper than those for either the coasting 
or the recollimation regions ((a CO ast) = 0.592 ± 0.007 and 
(^recoil) = 0.594 ±0.008, respectively). The mean for the de- 
celeration region, (adecei) = 0.618 =b 0.008, is in the middle 
of this range. While both the random and calibration errors 
contribute to the observed dispersion, any systematic error 
in the flux-density scale does not. We expect the rms scale 
error in spectral index to be < 0.02 (the value corresponding 
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Figure 6. Plots comparing spectral indices at different locations 
in the same jet. (a) High-emissivity and deceleration, (b) High- 
emissivity and coasting, (c) High-emissivity and recollimation. (d) 
Deceleration and coasting, (e) Deceleration and recollimation. (f) 
coasting and recollimation. Blue squares: main jets; red triangles: 
counter-jets. The error bars correspond to random errors a r only: 
calibration errors would be expected to move the point parallel 
to the line of equality, which is plotted for reference. 



to independent scale errors at 1.4 and 4.9 GHz as given by 
IPerlev Butler! I2013T ). 

More accurate measures of the change in spectral in- 
dex with distance are given by the differences for individ- 
ual jets, since any scale and calibration errors then cancel 
out. We plot the combinations of spectral indices at the 
four standard locations in Fig. [6j All of the combinations 
show spectral flattening with increasing distance except for 
the plot of a r ecoii against a CO ast (Fig- EF)- I n particular, all 
of the jets have spectra which flatten between the high- 
emissivity region and both the coasting and recollimation 
regions (Figs \6jp> and c). The mean differences are listed in 
Table [6] By this measure, the spectral flattening between 
high-emissivity and deceleration regions is not very signif- 
icant (2.5a) and that between coasting and recollimation 
regions is consistent with zero. The remaining differences 
show significant flattening (>4.5cr). The most significant dif- 
ference is between the high-emissivity and recollimation re- 
gions, for which (o^e — cWoii) = 0.067 ± 0.006, consistent 
with a constant spectral-index change for all of the 15 jets 
with measurements at these locations (Fig. [(J). The particu- 
larly good agreement between the integrated spectral indices 
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Figure 7. Comparison between the main and counter-jet spec- 
tral indices, ckj and cn c j for the four fiducial regions, (a) high- 
emissivity; (b) deceleration; (c) coasting; (d) recollimation. The 
error-bars represent the random errors, a T . 

measured for the outer two regions in the same jet, despite 
the evidence for slight local variations in spectral index asso- 
ciated with total-intensity enhancements (Section 14. 5p sug- 
gests that the internal errors cannot be much greater than 
we assume. 

The internal consistency of the measurements is ex- 
tremely good: all of the jets for which we have accurate 
spectral indices show the same trend. There are a few dis- 
crepant points but these all represent faint counter-jets for 
which the errors in our spectral index estimates are large, so 
the discrepancies are not significant (Fig. [6}. The rms spec- 
tral indices predicted by the combination of random and 
calibration errors, cr pre d, are also close to the observed val- 
ues (Table [5}. More formally, we find that our measurements 
and error model are consistent with all jets having the same 
spectral indices at the standard locations: this hypothesis is 
excluded only at the 73%, 40%, 11% and 74% confidence lev- 
els by a x 2 test for the high-emissivity, deceleration, coasting 
and recollimation regions, respectively. The validity of these 
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confidence levels is, however, almost entirely dependent on 
the accuracy of our estimate of calibration error: if we have 
overestimated this, then some intrinsic dispersion in spectral 
index would be required, albeit limited to <0.03 - 0.05 rms 
even for perfect data. 

Fig. [6] also shows strong correlations between the spec- 
tral indices measured in the same jet at different locations. 
The significance levels of these correlations, derived from 
the Spearman rank test, are given in Table [6] Four of the 
correlations are significant at >99.7% confidence; the excep- 
tions are those between a^e and a coas t or adecei- There are 
two possibilities: the correlations result from residual cal- 
ibration errors or the spectral-index evolution of the jets 
follows a roughly similar trend, but with source-dependent 
offsets. These cases are not mutually exclusive and are dif- 
ficult to distinguish, but we note that our assumed 3% cali- 
bration error corresponds to an rms spread in spectral index 
of ~0.03 between 4.9 and 1.4 GHz. This is very close to the 
rms spread along the line of equality in any of the panels of 
Fig. [6] If our assumed calibration errors are accurate, then 
we cannot exclude the possibility that all of the jets have 
identical spectral indices at the fiducial locations. 

4.4 Spectral differences between main and 
counter-jets 

We find no statistically significant differences between the 
spectral indices of the main and counter-jets at any of the 
four fiducial locations. The individual values are plotted in 
Fig. and the mean differences and their errors are given 
in Table [7] As noted earlier, the only source which shows 
a consistent offset between the two spectral indices is M 84 
(Fig. Uh), but this is almost certainly due to imaging arte- 
facts (Appendix [5j). The correlations between spectral in- 
dices for the main and counter-jets at the fiducial locations 
are not significant (<68% using the Spearman rank test; 
Table Q. 

If the spectrum is not a pure power law and the jets 
are relativistic, then systematic differences are expected be- 
tween observed spectral indices of approaching and receding 
jets. We return to this point in Section f5. 3 1 

4.5 Transverse and local spectral-index variations 

In detailed studi es of 3 C31, NGC315 and 3C296 
(Laing et al. 2006a .b, 2008b), we drew attention to two ad- 
ditional features of the spectral-index distributions. 

Firstly, the spectru m of the main jet i n NGC 315 flat- 
tens away from its axis ([Laing et al. 2006b) as well as with 
distance from the AGN; the same effect is visible in its 
counter-jet (Fig. [2J)) an d, at a low level, o n one side of 
each of the jets in 3C 31 (jLaing et alJl2008bl Fig. Et of the 
present paper). We have found no similar systematic pat- 
terns in the rest of the sample. It is doubtful, however, that 
we could have reliably detected such a subtle effect in sources 
with steeper-spectrum lobe emission surrounding the jets 
given the simplicity of our linear subtraction algorithm (Sec- 
tion 12. 3|) . Of the five sources where lobe emission is unde- 
tectable over the region of interest, NGC 315 and (at a much 
lower level) 3C 31 show the spectral flattening away from the 
axis mentioned earlier, 0326+39 and 3C 449 have no system- 
atic gradients and 1553+24 is not well enough resolved. The 



pronounced transverse spectral gradient in NGC 315 there- 
fore seems to be unusual, although we cannot rule out similar 
effects at a lower level in other sources. This gradient extends 
over many beam areas at high signal-to- noise (Fig. [2J)). In 
contrast, the narrow rims of steeper-spectrum emission seen 
close to the cores in NGC 193, 3C31, 3C270 and 3C 296 
(Figs c, g, i) occur in regions where the jets are nar- 
row and faint, and are likely to result from small calibration 
or deconvolution errors affecting the bright cores. Narrow 
spectral-index features at the faint edges of jets at larger 
distances from the nucleus are also suspect, particularly in 
cases where significant lobe emission has been subtracted 
(e.g. the counter-jet in 0206+35, Fig. Hi). 

Secondly, we identified "arcs" - narrow enhancements 
in total and polarized intensity crossing the jets - in 3C 31 
and 3C 296. The two brightest arcs in 3C31 have slightly 
Aa < 0.05) flatter spectra than the surrounding emission 
Laing et al.l 12008 bh . In 3C 296, the spectrum is flatter in 
both jets roughly where the more distinct arcs are seen, but 
the arcs are n ot recognisable indi vidually on the spectral- 
index images (jLaing et al.1 l2006ah . The remaining small- 
scale fluctuations visible in Fig. [2] are consistent with noise 
and calibration errors. 



5 DISCUSSION 

5.1 A common spectral-index profile? 

Our jet models suggest that the flaring region is an ap- 
proximately homologous structure, in the sense that the 
width and fiducial distances for emissivity and velocity all 
scale with the recollimation distance despite the wide range 
in physical size (the recollimation distance measured along 
the jet axis, ro, ranges from 1.5 kpc in M84 to 35kpc in 
NGC 315). We will discuss these relations in detail else- 
where (Laing & Bridle, in preparation), but the scaling is 
already evident from Table [3] The spectral indices at our 
standard locations in the flow show very little dispersion be- 
tween sources, so they must also follow this scaling, as can be 
illustrated directly by plotting the spectral index as a func- 
tion of distance normalized to the size of the flaring region, 
z/zq. In order to remove the confusing effects of scale and 
calibration errors, we evaluated the differences Aa between 
the binned spectral indices in Fig. [4] and a reference value in 
the same jet. We took the reference to be a C oast (cWoii for 
3C31 and 3C449, where the former is undefined), forcing 
the mean profile for each jet to have Aa = at the refer- 
ence location. The results are shown in Fig. [8] There is a 
clear trend in the combined data: on average, the spectral 
index falls from Aa w 0.06 at z/zo ~ 0.2 (as close to the 
flaring point as we can measure) to Aa « at z/zo « 0.8 
(the reference location). This confirms our new result: the 
variations of spectral index appear to scale with the size of 
the flaring region in the same way as variations of velocity 
and emissivity. 

Jets evolve in cross-section and velocity as a conse- 
quence of mass injection and prop agation in external pres - 
sure and density gradients (e.g. lLaing &; Bridle! l2Q02bh . 
Their emissivities are in turn governed by the balance of 
particle acceleration and energy-loss processes. The most 
natural way to ensure a common scaling for these quanti- 
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Figure 8. Plot of Aa against normalized distance from the AGN. 
Aa is the difference between the spectral index of a bin from Fig.@] 
and its mean value in the coasting region (or the recollimation 
region, if the former cannot be measured; see Table 0. Distance 
from the AGN, z, is normalized by the recollimation distance, 
zo- Blue squares: main jets; red triangles: counter-jets. Data at 
z/zo > 1.5 are available only for a few jets and are not plotted. 
0755+37 (which has no measured reference spectral index) is also 
omitted, as are data from upstream of the flaring point in 3C 270 
and 3C449. 
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Figure 9. Plot of Aa against jet flow velocity. Aa is the difference 
between the spectral index of a bin from Fig.0]and its mean value 
in the coasting region (or the recollimation region, if the former 
cannot be measured; see Table 0. The velocity /3o.5 is that for a 
streamline mid- way between the axis and the edge of the jet, from 
our model fits. Blue squares: main jets; red triangles: counter- 
jets. We have not plotted the inner two bins for the main jet of 
3C 270, as these include significant emission from upstream of the 
flaring point. 0755+37 (no reference spectral index) and 3C449 
(no velocity model) are also omitted. 



ties is for the acceleration mechanism to depend on flow 
speed, as we now discuss. 

5.2 A relation between spectral index and flow 
speed? 

Given that we model the jets as decelerating flows, and that 
the spectral index decreases monotonically with distance 
from the AGN in the regions studied here, there must be 
some relation between a and speed. In order to quantify 
this relation, we have chosen a representative speed /?o.5 for 
each source, determined from our model for a streamline 
half way between the axis and boundary of the jet. Spectral 
index difference (as in Section I5J"]) is plotted against /3o.5 in 
Fig. [9] There is a clear separation: all points with /3o.5 > 0.5 
have Aa > 0, whereas the mean A a at lower velocity is 
close to 0. 

An upper bound to the sound speed in the jet is that 
for an ultrarelativistic plasma, /3 S = 3 _1//2 = 0.58. /3 S may 
be significantly lower than this limit, particularly if the jet 
contains appreciable numbers of non-relativistic protons (en- 
trained or part of the initial composition), and may change 
system atically with position in the jet (|Laing &; Bridle! 
l2002bh . One possibility is then that the jet composition is 
dominated by ultrarelativistic leptons and magnetic field 
with a minor baryonic component and that the spectral 
index is significantly higher than its asymptotic value of 
(a) =0.59 wherever the average flow speed is even slightly 
supersonic. For the 10 sources we have modelled, the mean 
velocity at the flaring point is /3 — 0.81 with an rms of 
0.08. This in turn corresponds to a generalized Mach num- 



ber of M = (r/3)/(r s /3 s ) — 1.95 for an ultrarelativistic 
plasma, where T s = (1 - /3 S 2 )" 1/2 . For 3C31, we found from 
a conservation- law analysis that the flow throughout the jet 
base was transonic everywhere, with the Mach number de- 
er easiri^^c^A^^^ to_^M ^ 1, despite the increasing mass 
load (Lai ng Brid le 2002b). Any shocks in this flow regime 
must be weak. 

An alternative possibility is that the sound speeds are 
much lower than the ultrarelativistic limit, in which case 
larger Mach numbers are possible for the flow. 

5.3 Constraints on particle acceleration 

We first estimate some fiducial numbers for magnetic fields, 
for the radiating electron Lorentz factors and for energy 
loss timescales. The on-axis magnetic-field strengths derived 
from our emissivity models are ^2—9 nT at the flaring point 
and ~ 0.3 — 2nT at recollimatiorQ. The spectrum of syn- 
chrotron emission from a single electron is quite broad, but a 
characteristic frequency is v « 7 2 ^ g , where 7 is the electron 
Lorentz factor and z/ g = eB/27rm e is the gyro frequency. The 
corresponding Lorentz factor ranges for electrons radiating 
at 1.4 - 8.5 GHz are 2000 <7< 10000 at the flaring point 
and 5000 < 7 < 30000 at recollimation. The synchrotron loss 
timescale is 

2 We assume that: the sum of energy densities of relativist ic par- 
ticles and magnetic field is a minimum (close to the equipartition 
condition); there are no relativistic protons; the filling factor is 
unity and the energy spectrum is a power law between electron 
Lorentz factors of 100 and 10 6 . 
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E/(-dE/dt) = 



3ra e c/i 
2a T B 2 j 



which ranges from 6 x 10 5 to 9 x 10 7 yr for emission at 
8.5 GHz. 

A simple picture in which relativistic electrons are ac- 
celerated in the high-emissivity region and then passively 
advected along the jet, suffering radiative and adiabatic 
losses as they go, is not consistent with the spectral flat- 
tening we observe: radiative losses cause the spectrum to 
steepen at high energies. Even if the radiating particles are 
advected down the jet in a constant magnetic field, then 
the observed frequency spectrum would appear to steepen 
with distance. In fact, the magnetic field is likely to de- 
crease away from the AGN as a consequence of flux- freezing 
in a expanding flow. The consequent rapid decrease may 
be partially compensated by deceleration, shear or turbu- 
lent amplification, but (as mentioned earlier) the equipar- 
t it ion field strength still falls by a factor of ~5 - 10 be- 
tween the flaring point and recollimation. At a fixed fre- 
quency, we therefore observe higher-energy electrons and 
the steepening would be more pronounced. In any case, 
the minimum loss timescale for X-ray synchrotron emis- 
sion is ^100 yr for our assumed field strengths and a pho- 
ton energy of 1 keV. This corresponds to a flow distance 
of ~ 30//3pc (much smaller than the size of a typical jet 
base), so distributed reacceleration of radiating electrons 
throughout the regions we observe in FRI jets is unavoid- 
able. Independent evidence for ongoing acceleration comes 
from modelling of radio brightness and polarization evolu- 
tion along the jets: additional electrons are also required to 
offset the effects of adiabatic losses even at energies where 
the e ffects of synchrotron cooling are small (|Laing &; Bridle! 
2004). We must therefore attribute the change in character- 
istic radio spectral index along the jets to the underlying ac- 
celeration mechanism rather than to the combined effects of 
synchrotron and adiabatic losses. Although downward cur- 
vature at high frequencies is seen in the broad-band spectra 
of jet bases (lEvans et aTll2005l: iHardcastle et al.ll2Q0ll. liooBI: 
Lanz et al.1 1201 ll: IPerlman fc Wilson! 120051: IPerlman et al.1 



2010l : IWorrall et alJl20ld : lHarwood fe Hardcastlell2012h . the 
break frequencies are typically > 10 12 Hz, well above the 
range considered here. 

We have measured the spatial variation of the slope of 
the radiation spectrum over fixed ranges of observed fre- 
quency. In order to characterize the particle-acceleration 
mechanism, we need to derive the slope of the energy spec- 
trum over a fixed (or, at least, known) energy range. Unless 
the spectrum is a perfect power law, this presents several 
complications. 

(i) There is a known correction factor of 1 + z for the 
source redshift. 

(ii) We infer relativistic flow speeds for the jets, so the 
ratio of observed and emitted frequencies (^obs/^em = D, 
where D is the Doppler factor) may differ significantly from 
unity. The emitted frequencies for the approaching and re- 
ceding jets also differ systematically, but can be calculated 
from our jet models. 

(iii) As mentioned earlier, the magnetic field strength is 
likely to decrease away from the AGN along the jet, so emis- 
sion at a fixed frequency is generated by higher-energy elec- 
trons at larger distances. The evolution of field strength with 



distance is subject to uncertain and mutually contradictory 
assumptions such as equipartition with relativistic particle 
energy density or flux freezing. 

The correction for redshift is small for these nearby sources 
(Table[T]) and it is constant over both jets, so we can ignore it 
by comparison with the other two effects, which we consider 
in turn. 

The angles to the line of sight derived from our models 
range from = 25° to — 76° and the maximum velocities 
over the regions we measure spectral indices are in the range 
0.75 < f3 < 0.95. For const ant- velocity, antiparallel flows, 
the Doppler factors for the approaching and receding jets 



Dj = [r(l -/3cos6>)] _1 
Dd = ^(l + ^cos^)]- 1 



If the energy spectrum deviates from a pure power law, but is 
the same at a given distance from the AGN in both jets, then 
we would expect a systematic dependence of spectral index 
difference a c j — ctj on the ratio of Doppler factors Dj/D C y 
The maximum effect for our data should occur close to the 
flaring point, where the velocities are highest, and is in prin- 
ciple very significant for our sample: the maximum predicted 
ratio (for 1553+24) is DJ D C] ~ 13. Two problems make the 
effect hard to observe: (a) the largest Doppler factor ratios 
occur at small angles to the line of sight, so projection makes 
it difficult to resolve the flaring point from the AGN and (b) 
if D c j is small, the counter-jet emission is faint and errors 
on its spectral index are large. As a consequence, we have 
only been able to measure spectral-index differences over 
the range 1 < DJD^ < 3.5 for any of the fiducial regions, 
and those close to the AGN are less well determined. We see 
no systematic trend in plots of a c j — aj against DJD^ for 
either the high-emissivity or coasting regions. 

If all jets have the same electron energy spectrum 
which flattens with increasing energy, then there is a nat- 
ural explanation for the observed spectral gradients, since 
the magnetic-field strength must decrease away from the 
AGN. It seems implausible, however, that this (rather sub- 
tle) spectral flattening would occur in precisely the same 
observed frequency range for all of the jets in our sam- 
ple, despite the range of field strengths and Doppler fac- 
tors. The flattening cannot continue over a large energy 
range without violating constraints from broad-band obser- 
vations (e.g. IHardcastle et al. l l2002U2005h . Finally, it is un- 
clear what mechanism could cause the flattening. This pos- 
sibility should be tested using accurate measurements of jet 
spectra over a much larger frequency range, but we consider 
it to be very unlikely. 

The particle acceleration mechanism (s) at work in FRI 
jets must therefore satisfy the following constraints: 

(i) The electron energy spectrum can be modelled as a 
pure power law with spatially-varying mean index (p) = 
2(a) + 1 restricted to the range 2.18 < (p) < 2.32 (Table 0. 
The characteristic electron Lorentz factor range probed by 
our observations is 2000 < 7 < 30000, given the assumptions 
made earlier. 

(ii) The energy spectrum flattens with distance from the 
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AGN, from (p) — 2.32 in the high-emissivity region to (p) — 
2.18 in the coasting and recollimation regions. 

(hi) By the argument of Section [5721 (Fig. [9]), p may be 
determined by the characteristic bulk flow speed: a change 
in {p) from 2.32 to 2.18 being associated with a deceleration 
from p « 0.8 to ^< 0.5. 

(iv) The mechanism responsible for the steeper energy 
spectrum in the high-emissivity regions must be capable 
of accelerating electrons to Lorentz factors of 7 ~ 10 8 , as 
demonstrated by the detection of X-ray emission in many 
cases (Table[3] Fig. [3]). X-ray and/or optical emission is typi- 
cally seen from the high-emissivity region of the brighter (ap- 
proaching) jet; in some cases it extends slightly beyond the 
end of the high-emissivity region, but never beyond the end 
of deceleration. In the well-observed brighter jets of 3C 31 
and NGC 315, the ratio of X-r ay to radio emission decreases 
with distance from the AGN ([Laing et al. 2008a). This de- 
crease coincides with the spectral flattening we observe in 
the radio band and with deceleration. 

(v) The high-frequency spectrum at larger distances from 
the AGN (where we infer flatter energy spectra) is not 
known in most cases. The approaching jet of 3C 31 is de- 
tected at 8 p w avelength out to a distance of ~25 arcsec 
(|Lanz et alj l201lh . implying that acceleration of electrons 
to 7 « 10 5 — 10 6 (but not much higher) can occur after 
deceleration and recollimation. 

(vi) Acceleration is distributed, rather than being asso- 
ciated with rare, localized, prominent brightness enhance- 
ments. 



5.4 Shock acceleration 

First-order Fermi acceleration is thought to be the domi- 
nant particle acceleration mechanism at collisionless MHD 
shocks a nd there is an extensive lit erature about this pro- 
cess fsee lSummerlin Sz Barin3l2012l for a comprehensive list 
of references). It is also the only process for which there are 
comprehensive predictions of the energy index under differ- 
ent physical conditions, as we now discuss. 



5.4-1 Non-relativistic shocks 

Early work (e.g. [Bell 1 19781 : iBlandford fe Ostrikerl [l978h es- 
tablished that a power-law energy spectrum with index 
p = 2 (a = 0.5) is produced by strong, non-relativistic 
shocks in the test-particle approximation. This value is in- 
consistent with our results. Steeper spectra can in principle 
be produced in weaker shocks: p — (r + 2)/(r — 1) , where r is 
the velocity compression ratio, and 1 < r < 4. There is then 
no obvious rea son for p to h a ve a n arrow range, however, as 
pointed out bv lYoung et al.1 (|2Q05h . 



5.4-2 The ultrarelativistic limit 

For ultrarelativistic parallel shocks (i.e. those in which the 
magnetic field is parallel to the shock normal) in which parti- 
cles experience frequent small-angle scattering, a power-law 
energy spectrum with p = 2.23 (a = 0.615) is produced. 
This result has been found b oth analytically (|Kirk et al.l 
2000) and numerically (e.g. iBednarz '& Ostrowskil 1 1998: 
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lEllison fe Doubldl2004l : ISummerlin fe Barindl2012h . The co- 
incidence between the predicted spectral index and the 
value of 0.62 we fo und for three sources in our earlier work 
(jLaing et alJ l2008a) led us to wonder whether the ultrarel- 
ativistic limit might be relevant at the flaring point. This 
was puzzling, as the asymptotic value is approached only 
for shock Lorentz factors r > 10 and we infer r < 2 down- 
st ream of the flaring poin t. On the revised flux-density scale 
of IPerlev Butlerl (2013), however, the inferred spectra are 
slightly but significantly steeper, with (a) = 0.66 =b 0.01 in 
the high-emissivity region for the full sample. Given this dis- 
crepancy and the velocities ft < 0.9 that we infer downstream 
of the flaring point in FRI jets, it now seems even less likely 
that the ultrarelativistic limit is relevant to FRI jet bases, 
although a spectral index of 0.62 is only marginally incon- 
sistent with our result if the rms error due to uncertainties 
in the flux-density scale is as high as 0.02. Closer to the nu- 
cleus, we have no good e s timat es of the flow Lorentz factor, 
although Perl man et al ] (|201lh estimate r ~ 4 — 5 for the 
M 87 jet at knot HST-1. 

5.4-3 Mildly relativistic shocks 

The case of mildly relativistic shocks has been ad- 
dressed using Monte Carlo methods in the test- 
particle limit for oblique shocks ([Ellison &; Double] [2004; 
ISummerlin fe Baring] l2012h and in the non-linear case for 
parallel shocks ([Ellison fe Doublell2002h . These studies show 
that mildly relativistic shocks can generate a wide variety 
of power-law slopes for the energy spectrum, depending on 
the obliquity of the field in the shocks and the nature of the 
scatt ering as well as the shock speed ([Summerlin &; Baring! 
|2012| ) . Our observations require a narrow range of energy in- 
dices, with a slight dependence on mean flow (and therefore 
shock) velocity. The case which comes closest to meeting 
these requirements is that of a mildly relativistic velocity 
upstream of the shock front in the Bohm limit of frequent 
small- angle scatterings with a mean free path A clos e to the 
electron gyro-radius r g ([Summerlin &; Baring] 120121 Figs 7 
and 8) . The dependence of p on field obliquity is also weak 
under these circumstances. The case of an upstream flow 
velocity /3 U = 0.71, a velocity compression ratio of r = 3.02 
and a sonic Mach num ber A4 — 2.6 gives 2.35 < p < 2.5, 
depending on obliquity ([Summerlin &; Baring||2012L Fig. 8); 
for (3 U = 0.1 and r = 4 , the range is 2 < p < 2.04 
([Summerlin &; Baring] 12012, Fig. 7), close to the limit of 
p = 2 for a non-relativistic shock. We suggest that these 
two cases might bracket the range of physical conditions in 
FRI jet bases, the velocities in the first case being slightly 
too fast for the flaring point, while those in the second 
case are too slow for the flow after recollimation, which 
appears to remain significantly relativistic in most cases. It 
is plausible that the upstream fluid velocities in the shock 
rest frame are in the appropriate range, for example if the 
shocks are stationary or moving slowly in the observed 
frame, and that there is a range of field obliquities. 

5.5 What are the particle acceleration 
mechanisms? 

The first-order Fermi process at mildly relativistic shocks 
in the limit A « r g appears capable of generating energy 
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spectra with the correct indices and of accounting (at least 
qualitatively) for their variations along the jets. One possi- 
bility, therefore, is that this is the sole acceleration mecha- 
nism. A potential difficulty, however, is the distributed na- 
ture of the acceleration required to power the X-ray emis- 
sion, which cannot be restricted to a few localized shock 
sites: instead, a network of shocks extending over a sub- 
stantial distance along the jet appears to be required. Such 
shock networks could be identified with the complex, non- 
axisymmetric brightness structure in the high-emissivity re- 
gion (e.g. Fig. |TJi) or with the arcs ob served crossing the 
jets at larger distances in some sources (|Laing et al.l [2006a. 
l2008bl ) . Efficient deceleration of the jet by entrainment prob- 
ably requires a trans onic flow in which any shocks would b e 
very weak, however (|Bicknelllll994l ; lLaing &; Bridldl2002bl ). 

The most likely alternative to a pervasive shock net- 
work is that there are two distinct acceleration mechanisms 
operating in these jets. 



(i) The first mechanism (which could again be diffusive 
shock acceleration) dominates wherever the flow velocity 
/3 > 0.5 and in particular in the high-emissivity region. It 
can accelerate electrons up to Lorentz factors 7 « 10 7 — 10 8 , 
enabling them to emit X-ray synchrotron radiation, and has 
a characteristic energy index p = 2.32 over the energy range 
we sample. 

(ii) A second mechanism with a lower characteristic en- 
ergy index p = 2.18 becomes important when the flow ve- 
locity falls below j3 ~ 0.5, an d dominates after dec e lerati on. 
Gradual shear acceleration (jRieger fe Duffvl I2004L l200d ) is 
a plausible candidate, as our jet models show that there 
are transverse velocity gradients in these regions. We origi- 
nally favoured this idea for NGC 315 on the grounds that the 
spectrum flattens away from the jet axis pr ecisely where we 
infer velocity gradients ([Laing et al.ll2006bl ). but this effect 
seems not to be general f Section 14. 5)) . Second-order Fermi 
acceleration by plasma turbulence is a possible alternative 
mechanism. 



The existence of two different acceleration mechanisms 
in FRI jet bases has also been suggested based on ev- 
idence from high-resolution observations of X-ray, opti- 
cal and r adio emission from the nearest examples. Cen- 
taur us Aj^^^^^^^]^^;^^^^^^^]^^) and 
M87 dPerlman fe Wilson! 120051 : IPerlman et al.1 l2Qllh . For 
Cen A. lGoodger et al.l ([201ol ) suggest that stationary, X-ray- 
emitting knots are formed when lumps of dense material 
(e.g. molecular clouds or stars with high mass-loss rates) 
penetrate the jet, resulting in strong shocks with long-lived 
particle acceleration up to high Lorentz factors. Conversely, 
knots emitting only in the radio band (which often have high 
speeds) would be formed by weak shocks where no particles 
are accelerated to X-ray emitting energies. The proportion 
of radio emission from the two processes changes systemati- 
cally with distance from the AGN. It is not known whether 
there is a systematic difference in the radio spectra for the 
two types of knot, as would be expected if the difference is 
related to the spectral gradients we have discussed. 



6 SUMMARY 

The principal observational results from our imaging of jet 
base spectra for a sample of 11 FRI radio galaxies over the 
frequency range 1.4 - 8.5 GHz are as follows. 

(i) We find a narrow range of spectral indices, a, with the 
mean values at our fiducial locations ranging from 0.66 to 
0.59. The dispersion between sources is consistent with our 
error model. 

(ii) The radio spectral indices decrease (i.e., the spectra 
flatten) with distance from the nucleus on the scales we have 
studied. 

(hi) This spectral flattening occurs between the flaring 
point and the end of rapid deceleration as determined from 
our modelling of the jets as decelerating relativist ic flows; 
thereafter the spectral index is constant until after the jets 
recollimate. 

(iv) The spectral flattening is highly significant: the av- 
erage decrease in spectral index between the high-emissivity 
region immediately downstream of the flaring point and just 
after the jet recollimates is 0.067 =b 0.006 and all jets for 
which we have data show spectral flattening between these 
two locations. 

(v) The mean spectral index for the high-emissivity re- 
gion is (a) = 0.662 ± 0.011; coasting and recollimation re- 
gions have (a) = 0.592±0.007 and 0.594±0.008, respectively. 

(vi) There are no significant differences between the spec- 
tral indices of main and counter-jets at the same distance 
from the nucleus. 

(vii) All of the jets show essentially the same variation 
of spectral index with distance from the nucleus when the 
distances are normalized to the observed recollimation scale. 

(viii) The steeper spectra close to the jet flaring points 
appear to be associated with bulk flow speeds f3 > 0.5. 

(ix) The measured values of a correspond to indices of 
p = 2.18 — 2.32 for the electron energy spectrum over a 
Lorentz factor range of 2000 - 30000, assuming equipartition 
magnetic fields. 

The range of observed spectral indices and the inferred de- 
pendence on velocity could result from first-order Fermi ac- 
celeration by mildly relativistic shocks, in the limit that the 
scattering mean free path is close to the electron gyro-radius. 
Such a model would require a network of shocks in order to 
produce distributed particle acceleration. 

Alternatively, there may be two acceleration mecha- 
nisms. The first can accelerate electrons to Lorentz factors 
of 10 7 — 10 8 , is dominant for bulk flow speeds f3 > 0.5 close 
to the flaring point and has p = 2.32 at low energies. The 
second dominates at slower flow speeds and thus at larger 
distances, has p — 2.18 and a lower maximum energy. Shock 
acceleration is (again) a good candidate for the first mecha- 
nism 

In order to distinguish between these (and other) al- 
ternatives, accurate, well-sampled, spatially-resolved spec- 
tra over a very wide frequency range (ideally from low- 
frequency radio to X-ray) are required. More detailed pre- 
dictions of the energy spectra produced by various accelera- 
tion mechanisms, especially for Fermi acceleration by mag- 
netized, mildly relativistic shocks in the Bohm limit, would 
also be valuable. 



© 0000 RAS, MNRAS 000, 000-000 



FR I jet spectra 19 



ACKNOWLEDGEMENTS 

We thank the referee, Eric Perlman, for a careful reading of 
the paper. 



REFERENCES 

Abdo A.A., et al., 2011, ApJ, 727, 129 

Abreu P., et al., 2010, Astroparticle Physics, 34, 314 

Baars J.W.M., Genzel R., Pauliny-Toth U.K., Witzel, A., 

1977, A&A, 61, 99 
Bednarz J., Ostrowski M., 1998, Physical Review Letters, 

80, 3911 

Bell A.R., 1978, MNRAS, 182, 147 

Bicknell G.V., 1994, ApJ, 422, 542 

Birk G.T., Lesch H., 2000, ApJ, 530, L77 

Blandford R.D., Ostriker J.P., 1978, ApJ, 221, L29 

Bridle A.H., Baum S.A., Fomalont E.B., Parma P., Fanti 

R., Ekers R.D., 1991, A&A, 245, 371 
Canvin J.R., Laing R.A., 2004, MNRAS, 350, 1342 
Canvin J.R., Laing R.A., Bridle A.H., Cotton W.D., 2005, 

MNRAS, 363, 1223 
Cornwell T.J., Evans K.F., 1985, A&A, 143, 77 
Croston J.H., Birkinshaw M., Conway E., Davies R.L., 

2003, MNRAS, 339, 82 
de Vaucouleurs G., de Vaucouleurs A., Corwin H.G. Jr., 

Buta R.J., Paturel G., Fouque P., 1991, Third Reference 

Catalogue of Bright Galaxies, Springer, New York 
Ellison D.C., Double G.P., 2002, Astroparticle Physics, 18, 

213 

Ellison D.C., Double G.P., 2004, Astroparticle Physics, 22, 
323 

Evans DA., Hardcastle M.J., Croston J.H., Worrall D.M., 

Birkinshaw M., 2005, MNRAS, 359, 363 
Falco E.E., et al., 1999, PASP, 111, 438 
Fanaroff B.L., Riley J.M., 1974, MNRAS, 167, 31p 
Feretti L., Perley R., Giovannini, G., Andernach H., 1999, 

A&A, 341, 29 
Goodger J.L., et al., 2010, ApJ, 708, 675 
Greisen E.W., Spekkens K., van Moorsel G.A., 2009, AJ, 

137, 4718 

Guidetti D., Laing R.A., Murgia M., Govoni F., Gregorini 

L., Parma P., 2010, A&A, 514, A50 
Hardcastle M.J., Birkinshaw M., Worrall D.M., 2001, MN- 
RAS, 326, 1499 
Hardcastle M.J., Worrall D.M., Birkinshaw M., Laing, 

R.A., Bridle A.H., 2002, MNRAS, 334, 182 
Hardcastle M.J., Worrall D.M., Birkinshaw M., Laing R.A., 

Bridle A.H., 2005, MNRAS, 358, 843 
Hardcastle M.J., Worrall D.M., Kraft R.P., Forman W.R., 

Jones C, Murray S.S., 2003, ApJ, 593, 169 
Hardcastle, M. J., et al., 2007, ApJ, 670, L81 
Harris D.E., Finoguenov A., Bridle A.H., Hardcastle M.J., 

Laing R.A., 2002, ApJ, 580, 110 
Harwood J.J., Hardcastle M.J., 2012, MNRAS, 423, 1368 
Katz-Stone D.M., Rudnick L., 1997, ApJ, 488, 146 
Katz-Stone D.M., Rudnick L., Anderson M.C., 1993, ApJ, 

407, 549 

Kharb P., O'Dea CP., Tilak A., Baum S.A., Haynes E., 
Noel-Storr J., Fallon C, Christiansen K., 2012, ApJ, 754, 
1 



Kirk J.G., Guthman A.W., Gallant Y.A., Achterberg A., 

2000, ApJ, 542, 235 
Laing R.A., 2002, MNRAS, 329, 417 
Laing R.A., Bridle A.H., 2002a, MNRAS, 336, 328 
Laing R.A., Bridle A.H., 2002b, MNRAS, 336, 1161 
Laing R.A., Bridle A.H., 2004, MNRAS, 348, 1459 
Laing R.A., Bridle A.H., 2012, MNRAS, 424, 1149 
Laing R.A., Bridle A.H., Cotton W.D., Worrall D.M., 
Birkinshaw M., 2008a, Rector T.A., De Young D.S., eds, 
ASP Conf. Series 386, Extragalactic Jets: Theory and Ob- 
servation from Radio to Gamma Ray, ASP, San Francisco, 
p. 110 

Laing R.A., Bridle A.H., Parma P., Feretti L., Giovannini 

G. , Murgia M., Perley, R.A., 2008b, MNRAS, 386, 657 
Laing R.A., Canvin J.R, Bridle A.H. Hardcastle, M.J., 

2006a, MNRAS, 372, 510 
Laing R.A., Canvin J.R, Cotton W.D., Bridle A.H., 2006b, 

MNRAS, 368, 48 
Laing R.A., Guidetti D., Bridle A.H., Parma P., Bondi M., 

2011, MNRAS, 417, 2789 
Lanz L., Bliss A., Kraft R.P., Birkinshaw M., Lai D.V., 

Forman W.R., Jones C, Worrall D.M., 2011, ApJ, 731, 

52 

Leahy J.P., Perley R.A., 1991, AJ, 102, 537 
Lemoine M., Pelletier G., 2003, ApJ, 589, L73 
Miller N.A., Ledlow M.J., Owen F.N., Hill J.M., 2002, AJ, 
123, 3018 

Morganti R., Fanti C, Fanti R., Parma P., de Ruiter H.R., 

1987, A&A, 183, 203 
Morganti R., Parma P., Capetti A., Fanti R., de Ruiter 

H. R., Prandoni I., 1997, A&AS, 126, 335 

Ogando R.L.C., Maia M.A.G., Pellegrini P.S., da Costa 

L.N., 2008, AJ, 135, 2424 
Parma P., de Ruiter H.R., Capetti A., Fanti R., Morganti 

R., Bondi M., Laing R.A., Canvin J.R., 2003, A&A, 397, 

127 

Perley R.A., Butler B.J., 2013, ApJS, in press, 
larXiv:1211. 13001 

Perlman E.S., Biretta J.A., Zhou F., Sparks, W.B., Mac- 

chetto F.D., 1999, AJ, 117, 2185 
Perlman E.S., Wilson, A.S., 2005, ApJ, 627, 140 
Perlman E.S. et al., 2010, ApJ, 708, 171 
Perlman E.S. et al., 2011, ApJ, 743, 119 
Petrosian V., 2012, Space Science Reviews, 173, 535 
Rieger F.M., Duffy P., 2004, ApJ, 617, 155 
Rieger F.M., Duffy P., 2006, ApJ, 652, 1044 
Smith, R.J., Lucey, J.R., Hudson, M.J., Schlegel, D.J., 

Davies, R.L., 2000, MNRAS, 313, 469 
Stawarz L., Ostrowski M., 2002, ApJ, 578, 763 
Summerlin E.J., Baring M.G., 2012, ApJ, 745, 63 
Trager S.C., Faber S.M., Worthey G., Gonzalez J.J., 2000, 

AJ, 119, 1645 
Weiland J.L., et al., 2011, ApJS, 192, 19 
Worrall D.M., Birkinshaw M., Hardcastle M.J., 2001, MN- 
RAS, 326, L7 

Worrall D.M., Birkinshaw M., Laing R.A., Cotton W.D., 

Bridle A.H., 2007, MNRAS, 380, 2 
Worrall D.M., et al., 2008, ApJ, 673, L135 
Worrall D.M., Birkinshaw M., O'Sullivan E., Zezas A., 

Wolter A., Trinchieri, G., Fabbiano G., 2010, MNRAS, 

408, 701 

Young A., Rudnick L., Katz D., Delaney T., Kassim N.E., 



© 0000 RAS, MNRAS 000, 000-000 



20 R.A. Laing et al 

Makishima K., 2005, ApJ, 626, 748 



APPENDIX A: ARCHIVAL VLA 
OBSERVATIONS OF 0326+39 AND 1553+24 

The observations and data reduction for the images of 
0326+39 and 1553+24 in the 4.5 - 5.0 and 1.3 - 1.7 GHz 
bands have not been documented elsewhere, so we sum- 
marize them here. Data were retrieved from the VLA 
archive, calibrated and imaged using the MPS package, 
following st a ndard procedures essentially as described by 
lLaing et al.l (|201lh . 3C48 or 3C 286 were used as primary 
amplitude calibrators. All sources were self-calibrated and 
data from multiple array configurations were combined. The 
final images were m ade using a multi-scale clean algorithm 
(jGreisen et al.ll2QQ9h . 

For 1553+24, we chose not to make images at 
1464.9 MHz, as the available A-configuration data at this 
frequency were badly affected by interference. We also anal- 
ysed other archival datasets at 1665 and 1415 MHz, but these 
were not observed in all relevant configurations and did not 
help to improve the accuracy of the final spectral-index im- 
age. 



APPENDIX B: 
SOURCES 



NOTES ON INDIVIDUAL 



NGC 315 A detail ed study of the sp ectral-index distri- 
bution was given bv lLaing et al.l ([200 6b ). The spectrum of 
the jets is unusually flat, although less extreme than in our 
earlier work after correction of the flux-density scale. 

3C 31 A detaile d study of t h e spect ral index distribution 
was presented bv lLaing et al.l (|2008bl ). The coasting region 
is not defined, as the deceleration region extends slightly 
beyond recollimation. 

0206+3 5 Spectral-index im ages of the whole source were 
shown by lLaing et alJ (| 20 111 ) and for the jets alone (with 
a less stringen t blan king criterion than used here) by 
lLaing fe Bridlel (|2Q12h . A residual calibration error caused 
a high-frequency ripple in spectral index transverse to the 
jet axis, visible in the brighter parts of the jets (Fig. [3jl). 
The amplitude of this ripple is roughly Act ~ 0.01, and we 
do not believe that it has any systematic effect on our re- 
sults. The high-emissivity region is not well enough resolved 
from the core for us to determine accurate spectral indices. 
We have included the part of the jet and counter- jet em is- 
sion that we model as backflow ( Laing &; Bridlel 120121 ) m 
the binned spectral-index profiles and averages, since these 
are very hard to separate from the outflow component at 
1.2-arcsec r esolution. 

0326+39 iBridle et all (|l99lh show spectral-index profiles 
derived from lower-resolution images at 1.5 and 4.9 GHz. 

0755+37 Spectral-index images of the whole source and 
for the jets have been presented elsewh ere, as for 0206+35 
(jLaing et al.ll201ll : lLaing fc Bridlel 1 20121 ). Only the main jet 
has enough unblanked emission for us to derive a reliable 
spectral index. We considered using lower-resolution images 
to improve the signal-to-noise ration, but lobe subtraction 
was problematic. Althou gh we model the jets as a mixture 
of outflow and backflow ( Laing &; Br idle 2012), the spectral 



Table Al. Journal of VLA observations. (1) source name; (2) 
date of observation; (3) VLA configuration (H denotes a hybrid 
configuration, as described bv lBridle et"aH ri991^: (4) and (5) fre- 
quencies (MHz); (6) bandwidth (MHz); (7) duration (min); (8) 
key to proposal code and reference. 
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Proposal codes and r eferences : 

1 BRID dBridle et al.lll99ll ) 

2 AR386 ( unpublished) 

3 AM154 dMorganti et al.lll987l ) 

4 AM222 ( unpublished) 

5 AM364 (iMorganti et al.ll 19971 ) 

6 AR402 dYoung et al.ll2005h 



indices in the present paper are dominated by the outflow 
component. 

3C 270 A detailed discussion of the radio structure, spec- 
tral index and magnetic-field structure of 3C 270 will be 
given by Laing, Guidetti & Bridle (in preparation). The av- 
erage spectral index of the jets of 3C 270 is unusually high. 
This is the only source where we can measure the spectral 
index accurately upstream of the flaring point. The abrupt 
steepening of the spectrum of the main jet from a « 0.55 
to a « 0.7 between distances of ^9 and ^12arcsec from the 
nucleus is seen most clearly in the unbinned profile (Fig. [3^). 
Given that our model fits give an unusually large range for 
the flaring-point position (6.5 - 11.5arcsec in projection), 
it is likely that the steepening occurs exactly at the flaring 
point. The location of the spectral steepening is blurred by 
the binning in Fig. Hfg). 

M 84 There are significant artefacts on the 1413-MHz 
image of M 84, caused by a combination of poor spatial- 
frequency coverage, the proximity of the bright confusing 
source M87 and the displa cement of the po inting centre in 
some of the observations ([Laing et al.ll201ll ). There are rip- 
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pies of amplitude Act ^0.1 and wavelength r anging from 
^3 t o ^15arcsec on the spectral-index image ([Laing et al.l 
l201ll . fig. 9c). These are aligned in the NS direction (roughly 
parallel to the jets). 

3C 296 We have given a comprehensive disc ussion of the 
spectral-ind ex distribution of 3C 296 elsewhere ([Laing et al.1 
l2006aLl201l| y. 



1553+24 fYoung et al.1 (|2005h showed images of this 



source (we include some of their visibility data) but did not 
discuss its spectrum quantitatively. The high-emissivity re- 
gion is too close to the core for accurate spectral-index mea- 
surements to be made; this is the one source where it does 
not overlap with the deceleration region. The signal-to-noise 
ratio in the counter-jet is low. 

3C 449 The spectra o f the j ets in this source were mea- 
sured by iFeretti et al.l (119991 ) and studied in detail by 
iKatz-Stone &; Rudnickl (|l997h . Our measured spectral in- 
dices are consistent with the mean value (a) = 0.58 =b 0.02 
betw een 4.985 and 8.4 G Hz with 2.5-arcsec resolution quoted 
by IFeretti et al.1 ^999) from the same visibility data once 
we account for recalibration. Note that the steep-spectrum 
emission surrounding the jet s does not start until well out - 
side our region of interest (Ka tz-Stone &; Rudnickl [l997h . 
The inner jets are very symmetrical, indicating that they 
are almost precisely in the plane of the sky. For this reason, 
we cannot model the jet velocity field, and have no good 
estimate of the deceleration scale, although by analogy with 
other sources it is likely to be < 10 arcsec. 



© 0000 RAS, MNRAS 000, 000-000 



